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Abstract 

Starting from the Bethe Ansatz solution of tlie open integrable spin-1 XXZ quan- 
tum spin chain with diagonal boundary terms, we derive a set of nonlinear integral 
equations (NLIEs), which we propose to describe the boundary supersymmetric sine- 
Gordon model BSSG^ with Dirichlet boundary conditions on a finite interval. We 
compute the corresponding boundary S matrix, and find that it coincides with the 
one proposed by Bajnok, Palla and Takacs for the Dirichlet BSSG"*" model. We de- 
rive a relation between the (UV) parameters in the boundary conditions and the (IR) 
parameters in the boundary S matrix. By computing the boundary vacuum energy, 
we determine a previously unknown parameter in the scattering theory. We solve the 
NLIEs numerically for intermediate values of the interval length, and find agreement 
with our analytical result for the effective central charge in the UV limit and with 
boundary conformal perturbation theory. 
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1 Introduction 



Much can be learned from studying quantum field theories in finite volume. This is particu- 
larly true for 1 + 1-dimensional integrable QFTs, for which there are effective descriptions in 
both the infrared (infinite size) and ultraviolet (zero size) limits, namely, massive factorized 
scattering theory [1, 2, 3] and conformal field theory (CFT) [4, 5], respectively. Moreover, at 
least for some examples, there also exist effective descriptions in terms of certain nonlinear 
integral equations (NLIEs) for general system size. 

A case in point is the sine-Gordon (SG) model on a circle (periodic boundary conditions) 
[6]- [8], and on an interval with either Dirichlet [9, 10] or general integrable [11] boundary 
conditions. NLIEs have been obtained in these papers from Bethe Ansatz solutions of corre- 
sponding critical spin-| XXZ chains [12]-[15], with a mass scale introduced by means of an 
alternating inhomogcncity parameter (±A). Among the quantities that have been computed 
from these NLIEs arc bulk and boundary S matrices (IR limit), bulk and boundary energies, 
and central charge and conformal dimensions (UV limit). Moreover, Casimir energies in the 
near UV region obtained by numerically solving the NLIEs agree well with those obtained 
using the truncated conformal space approach (TCSA) [16, 17] and conformal perturbation 
theory. 

In the NLIE approach one generally deals with three sets of parameters: the UV pa- 
rameters appearing in the action, the IR parameters appearing in the S matrices, and the 
lattice parameters in terms of which the NLIE is initially formulated. (See Figure 1.) In 
principle, by matching the UV and IR limits of the NLIE with corresponding known results, 
one can deduce the "lattice UV" and "lattice ^ IR" relations, respectively. If the "UV 
^ IR" relation is also known, then the consistency of these three sets of relations can be 
checked. For the sine-Gordon model, for which the "UV ^ IR" relation for the boundary 
parameters has been determined [18], this consistency has been established for both the bulk 
and boundary parameters. 

Much less is known about the supersymmetric sine-Gordon (SSG) model [19]-[22], ^ 

£ = ^d.ip d^ip + ^t^rd,^ - ^ cos(/?(^)f'* + ^ cos2(/?(^) , (1.1) 

where is a real scalar field, \E' is a Majorana Fermion field, mo is a mass parameter, and 
e (0 ,47r) is a dimensionless coupling constant. Indeed, for the case of periodic boundary 
conditions, an NLIE was proposed only recently [23] , and derived in [24, 25] from the Bethe 
Ansatz solution of the integrable spin-1 XXZ chain [26, 27]. While the ground state of 

iWetake7°= ( ^, ] , 7I = ( ^ * ) , * = [ ^ ), with V and i/i real. 
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Figure 1: Three sets of parameters and their relations 



the critical spin- 



I- 1 chain is described by a sea of real Bethe roots, the ground state of the 



critical spin-1 chain is described by a sea of approximate "two-strings", i.e., certain complex 
conjugate pairs of Bethe roots. As a result, the familiar method [7] of deriving the NLIE, 
based on the Bethe Ansatz equations and the corresponding counting function, does not 
seem to work for the spin-1 case. Nevertheless, an NLIE can be derived [24, 25] from the 
model's T — Q equations, in a manner similar to the original approach [6] . Very recently, the 
periodic spin-1 XXZ/SSG NLIE for excited states was shown [28] to have the correct UV 
and IR limits. In particular, the bulk soliton S matrix [22] was obtained from the IR limit 
of the NLIE. 

Our goal has been to find an NLIE for the boundary supersymmetric sine-Gordon model 
on an interval with general integrable boundary conditions [29]- [32]. This field theory is of 
interest as a toy model with boundary that is both integrable and supersymmetric; and it 
may also have applications to superstring theory [33] . Since a Bethe Ansatz solution of the 
corresponding open spin-1 XXZ chain with general integrable boundary terms [34] is still 
not known, we focus here on the special case of diagonal boundary terms, whose solution is 
already available [35]. The homogeneous chain has the local Hamiltonian 



where the bulk terms Hn,n+i are those of the Fateev-Zamolodchikov [26] spin chain. 



7V-1 




(1.2) 
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n n n+1 ' 
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and Sn are spin-l generators of SU{2); and the diagonal boundary terms are give by 
b.t. = ^ sin(27) | - [cot 7y_ + cot(7y_ - 7)] + [cot ?7_ - cot(77_ - 7)] {S^f 

- [cot77+ + cot(r/+-7)]5^ + [cotr/+-cot(r;+-7)](5^)2}. (1.5) 

The bulk and boundary parameters are 7 and r)±, respectively. 

Wc propose that, in analogy with the spin-i XXZ/SG model, this open spin chain cor- 
responds to the boundary SSG model on an interval [a;_ , x+] with Dirichlet boundary con- 
ditions, 

where ijj and ip are the spinor components of the Majorana field These boundary condi- 
tions follow from the boundary action [30] in the limit that the boundary mass parameters 
tend to infinity. That the diagonal boundary terms (1.5) of the spin-l XXZ chain correspond 
to the Dirichlet boundary conditions (1.6) of the SSG model is consistent with the fact that 
and topological charge are conserved in the two models, respectively, and also that both 
models are integrable. 

There are actually two known sets of integrable supersymmetric boundary conditions for 
the boundary SSG model [30]. Following [32], we shall refer to the set (1.6) as Dirichlet 
BSSG^, and to the set with ip + ip = at both ends as Dirichlet BSSG^. 

We derive an NLIE for the Dirichlet BSSG+ model, circumventing (as in [25]) the difficul- 
ties posed by the ground-state sea of two-strings by identifying suitable auxiliary functions 
from the model's T — Q equations [35], and exploiting their analytic properties. By analyzing 
the IR limit of this NLIE, we compute the soliton boundary S matrix, which coincides with 
the one proposed for the Dirichlet BSSG^ model by Bajnok et al. [32]. We propose the 
"UV ^ IR" relation for the boundary parameters, for the special case of Dirichlet boundary 
conditions, on the basis of our NLIE and its UV and IR limits. By computing the boundary 
vacuum energy, we determine a previously unknown parameter in the scattering theory [32] . 
We solve the NLIEs numerically for intermediate values of the volume, and find agreement 
with our analytical result for the effective central charge in the UV limit and with boundary 
conformal perturbation theory, and confirm the UV-IR relation. 

The outline of this paper is as follows. In Section 2 we rederive the NLIE for the spin- 
I XXZ/sine-Gordon model with Dirichlet boundary conditions [9]. However, we use the 
method which we apply to the spin-l case (which differs from the approach used in [9]), 
and therefore, this serves as a vahiable warm-up exercise for the latter problem. In Section 
3 we turn to our main interest, the spin-l XXZ/supersymmetric sine-Gordon model with 
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Dirichlet boundary conditions. We sketch the derivation of the NLIE, relegating some of 
the details to Appendix B. In particular, we determine the boundary terms Pbdryid) and 
Py{0) which encode boundary effects. Section 4 is devoted to an analysis of the IR limit of 
this NLIE. In the course of computing the corresponding boundary S matrix, we determine 
the "lattice ■*-> IR" relation for the boundary parameters. In Section 5, we analyze the UV 
limit of our NLIE, and compare with the expected CFT result. In this way, we obtain a 
boundary "UV ^ lattice" relation, and therefore finally the "UV IR" relation for the 
boundary parameters. In Section 6 we compute the effective central charge of the Dirichlet 
SSG model using first-order boundary conformal perturbation theory, and compare with 
numerical NLIE results. We conclude in Section 7 with a discussion of our results and with 
some comments on various open problems. Some important technical details are explained 
in the appendices. 



2 Spin-^ XXZ/SG with Dirichlet boundary conditions 

We rederive here the NLIE for the spin-| XXZ/sine-Gordon model with Dirichlet boundary 
conditions [9]. However, in contrast to the familiar approach [7] used in [9], we do not 
introduce the counting function. Instead, we identify suitable auxiliary functions from the 
model's T — Q equation, and exploit their analyticity properties. We shall use the same 
method to treat the spin-1 XXZ/supersymmetric sine-Gordon model in Section 3. 



2.1 T — Q equation 

The transfer-matrix eigenvalues T{x) of the inhomogeneous open spin-| XXZ chain with 
Dirichlet boundary conditions satisfy the T — Q equation [14] 

T{x) = T^+\x) + T^-\x) , rW(x) = sinh(2x ± ij) B^^\x) (t>{x ± ^ ^2.1) 

2 Q[x) 

where 

(j){x) = sinh^(x- A)sinh^(x + A), 

= sinh(a; ± sinh(a; ± , 

M 

Q{^) — sinh(a; — v^) sinh(a; -|- v^) ■ (2.2) 

k=\ 

We denote the bulk parameter by 7 (0 < 7 < vr), and the boundary parameters by H±. 
Moreover, A is the inhomogeneity parameter which provides a mass scale; N is the number 
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of spins; and the M zeros Vk of Q{x) are the Bethe roots. Note that T{—x) — —T{x). 

For the homogeneous case A = 0, a local Hamiltonian is obtained from the first derivative 
of the transfer matrix T{x) [14]. However, for the inhomogeneous case A 7^ which we 
consider here, the definition of energy is less clear. We shall follow the prescription of 
Reshetikhin and Saleur [36], which implies 



E^-l\^\nT{x) 
a I ax 



x=A+n. 



In T(a;) 

ax 



x=A- 



where a is the lattice spacing, and g is given by 

One can verify that this definition has the correct A — > limit. 



(2.3) 



(2.4) 



2.2 Derivation of NLIE 



We define the auxiliary functions a{x) and a{x) by 



a{x) 



sinh(2a; + ^7) B^+\x) (p^x + ^) Q{x - i-f) 



sinh(2a; — ^7) \x) 4>{x — y) Q{x + ^7) ' 
The transfer-matrix eigenvalues can then be written as 



a{x) — a{—x) 



a{x) 



T{x) = smh{2x-i-^)B^-\x)4){x-^)^^^^^^A{x) 
= sinh(2a; + i7)sW(x)0(x+^)^^|^^^(x), 



where 



A{x) = 1 + a{x) , A{x) = 1 + a{x) 
The Bethe Ansatz equations are given by 

A{vk) = Q, k = l,...,M. 



(2.5) 



(2.6) 



(2.7) 



(2.8) 



We consider the ground state. For simplicity, we restrict the boundary parameters H± 
to the interval 



0<H±< 



TT 



7 



(2.9) 



5 



We argue in Appendix A.l that the boundary parameters should be further restricted to the 
range 



l-2<H^ + H.<^-^-2 
7 7 



(2.10) 



in order for the ground state to have M = N/2 real Bcthc roots and no holes, except 
for one hole at the origin. That is, T{x) does not have zeros near the real axis except for a 
simple zero at the origin. To remove this root, we define 



f{x) 



(2.11) 



where ^{x) is any function whose only real root is a simple zero at the origin, in particular 
/i(0) = , //'(O) 7^ 0, so that T{x) is analj^ic and nonzero (ANZ) when x is near the real 
axis. (We use a prime {') denote differentiation with respect io x.) It is convenient to 
introduce the compact notation 

= t-{x) %^A(x) = U{x) ^^A{x) , (2.12) 



Q{x) 



Q(x) 



where 



, , sinh(2a; ± i7) ^ , , ,, ij . 



Since \nT{x) is analytic near the real axis, Cauchy's theorem gives ^ 

0^ ^ dx [\nf{x)]"e 



(2.13) 



(2.14) 



where we choose the contour C as in Figure 2. 

Dividing the contour C = Ci + C2, where Ci and C2 have imaginary parts ±ie with e 
small and positive respectively, this integral can be written as 



^ [ dx [lnt+(x)]"e^*^^+ / 
Jci Jc 

+ [ dx []lit-{x)]" 



dx < In 




e''=^+ / dx [\nA{x)] e 
Jci 



" Akx 



6^^=^+ / dx {\n 
C2 



In terms of Fourier transforms defined along C2 and Ci by 



e*'^"+ / dx [\nA{x)\'e 
C2 



" Akx 



Lf"{k)= dx [Infix)]" e'^^ Cf"{k) = [ln/(a;)] , 

' C2 Jci 



(2.15) 



^Since T{x) grows exponentially with x for large x (as follows from (2.1)), [lnT(a;)]" for x 00. 
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\e 








-ie 




^2 



Figure 2: Integration contour 



respectively, we can rewrite 



dx < In 



Ci 



dx < In 



C2 




ikx 



ikx 



-LQ"{k) {e-^^ -e-^^) , 



LQ"{k) [e^'<--^^ - 1] , 

where we have used the periodicity 

Q{x) = Q{x — m), X G Ci, and (5(a; + 27) = + ^7 — ivr), x G C2 . 
Defining 



C{k) 



dx[\Yit+{x)\ e'"'' + / dx[\nt^{x)\ e 

Cl JC2 



" Akx 



we obtain 



LQ"{k) 



Trfc 

e 2 



4cosh(f )sinh((7r-7)|) 

This is the main consequence of analyticity. 

It follows from the definition of a{x) (2.5) that 

' Q{x - i-f) 
_Q{x + i'y) 



LA"{k) + CA"ik) + Cik) . (2.19) 



La"{k) = I dx Un 



e*fe^ + D{k) 



where D{k) is defined by 



D{k) = I dx Un 



sinh(2a; + 27) 5(+)(x) 0(x + |) 
sinh(2x - n) B(-){x) 0(a; - f ) 



Akx 



(2.16) 



(2.17) 



(2.18) 



(2.20) 



(2.21) 
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Using again the periodicity of Q{x), we obtain 

L^'{k) = LQ"{k) [e-^^ - e^T-")'^] + D{k) . 
Inserting (2.19) into (2.22) yields 

L^'{k) = G{k) \LA''{k) + £A"{k) \ + Crik) , 

where 

G(k) = sinh((vr-27)f) 
^ ^ 2cosh(f )sinh((7r-7)|) ' 

Crik) = G{k)C{k)+D{k). 



(2.22) 
(2.23) 

(2.24) 
(2.25) 



The evaluation of CT{k) is tedious but straightforward. To this end, we make use of the 
identities 



dx 



^ [In sinh(x-iQ;)]" 6*^=^ = 6 



" ikx — fe(a— riTT 



C2 



2ti 



k 



l-e 



-Trfe 



(2.26) 



where n is an integer such that < ^e{a — mr) < tt, and 

k 

— [lnsinh(2a;)J" e'""" = — 



/ — [lnsinh(2a;)]"e'*^^ = - 
Jc2 27r 1 



(2.27) 



e 2 



The result is (B.6) 
Crik) = 



^^^J n cos(Ak) sinh {{^H+ - 7r)|) + sinh {{-fH^ - 7r)|) 
{ cosh(f ) ^ 2cosh(f)sinh((7-7r)|) 

cosh(^)sinh((27-7r)|) 



cosh(f )sinh((7-7r)|) 



(2.28) 



The result (2.23) is the NLIE for the lattice sine-Gordon model with Dirichlet boundary 
conditions in Fourier space. Passing to coordinate space, and integrating twice, we obtain 



\na{x) 



- i2A^tan 



dx' G{x - + ie) \n.A{x' - ie) - / dx' G{x - x' - ie) \nA{x' + ie) 
( sinh — 

1 I 7 



\ cosh 

\ 7 



^ + i Pbdry{x) + in , 



(2.29) 



where G{x) is the Fourier transform of G{k) (2.24) 

1 f°° ^ 
G{x) = — dk 6-^^=^ G{k) , 
27r ./_^ 



(2.30) 



Pbdry{x) is given by 

Pbdryix) = dx'R{x')^-l dx'R{x'), (2.31) 
Jo 2 J_j. 

and R{x) is the Fourier transform of R{k), 

. ^ ^_^[ sinh ((7^+ - tt)!) + sinh {{^H^ - 7r)|) 
^ ^ 2cosh(|;)sinh((7-7r)|) 



cosh(^)sinh((27-7r)|) 



(2.32) 



cosh(f )sinh((7-7r)|) J ' 

The integration constant is explained in Section C.l. 

The continuum hmit consists of taking A ^ oo (leading to a simphfication in the driving 

_i / sinh 2^ \ _eA 

term, — i2A'"tan I A j ~ —lANe t sinh y), together with — > oo and lattice spacing 
a — > 0, such that the interval length L and the soliton mass m are given by 

L^Na, m^-e'^f. (2.33) 
a 

The driving term therefore becomes —i2mLsinh9, where the renormalized rapidity 6 is 
defined as 

e^—. (2.34) 

7 

The resulting NLIE 

de' G{e-e' + is) in ^{e' -is) - / do' G{e-e' - ie) in A{e' + is) 

-OO J —OO 

— i2mL sinh 9 + i Pbdryi^) + ^tt , (2.35) 
where we have defined 

£ = a(^) = a(^), A{9) = A{^), P,,^y{9) = P,,ry{-) , G(^) = ^G(^) ,(2.36) 

7 TT TT TT TT TT 

agrees with [9] . 

2.3 Vacuum and Casimir energies 

Since our approach avoids introducing the counting function and the density of Bethe roots, 
the energy computation also differs from that of the conventional approach [7]. The main 
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idea is to express the energy in terms of LQ{k), and then make use of the consequence (2.19) 
of analyticity. 

The energy is given by (2.3). From the T — Q equation (2.1) and the fact 0(A) = 0, we 
see that 



d_ 
dx 



hiTix) 



dx 



\nT^^\x) 



x=A±n. 



Hence, the energy (2.3) can be written as 



9 d 



a dx 



E = _^_L<jinTW(a; + ^)-lnT^-^(a;-^) 



aj 27r 

where we used the fact (see (2.15)) 



X eC2. 



Since 



T^^\x ± ^) = sinh(2x ± 2i7) 



(f){x ± i'j) 



Q{xTf) 



we can compute the Fourier transforms e^^LT^^y{k), ^ 



e-^LT(-y{k) = e-'^LB(-y{k) + e-^''L(i)'{k) + e 
- 26-"^ sinh((7r-7)^)Lg'(A;), 



i-ik) 



e'^LTM'ik) 



{-ik) 

k k - 

+ 2e-'^sinh((7r-7)-)Lg'(A;). 



We can ehminate LQ<{k) using (2.19) and the fact Lf<{k) = ^Lf"{k), 
1 



LQ'{k) 



irk 

e 2 



LA"{k)+ CA"{k) + C{k) 



(2.37) 



(2.38) 

(2.39) 
(2.40) 



(2.41) 



{-ik) 4cosh(^) sinh ((tt - 7)|) 

Trfc 

, ' 1 + ^'i^) + LB^Hk) - IBM' 

4cosh(^) smh ((tt - 7)|) 



{k) 



+ 



L(l)'{k) + [e 



ik ( J TT \ 7, 



2TrHk)_ 
i-zk) 



-27ri}, (2.42) 



^We assume here that < 7 < 7r/2. 
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where we have passed to the second hne using (see (B.2), (2.21)) 



c{k) = LB(-y'{k)-LB(+y'{k) + 



' _ 7*1 ( X 

e 2 — eU 



L(P"{k) 



+ 



e 2 -eV2 2;'= 



(2.43) 



We obtain 



e^UrM'{k)-e-'-^lJT^y{k) 



1 



cosh 2^ 



+ tanh(^) 



2e-^ cosh((7 - ^)k)I4'{k) 



LA'{k) + CA'{k) - 2m 
27rcosh((7 - j)k) 



i sinh ^ 



+e-'^LBi^'{k) + e^LB^y{k) 



This is essentially the integrand in the expression (2.38) for the energy. 

Let us consider this result term by term. The L(f)'{k) is obtained as before (B.5) 



(2.44) 



-lAk-' 



m 
i-tk) 



Inserting the corresponding order contribution from (2.44) into (2.38), we obtain the bulk 
energy 



Ng_ 

ia 



dk (1 + e-^'^'') 



sinhf cosh((7- |)A;) 



Trk 



cosh ^ sinh ^ 



(2.45) 



This quantity is divergent in the continuum limit. We adopt the renormalization procedure 
[7] of discarding divergent terms and keeping only the (finite) terms that can be expressed 
in terms of the physical mass m given in (2.33). To this end, we discard the A-independent 
term, and evaluate the remaining integral by closing the contour in the lower half plane and 
selecting only the contribution from the pole at k = We arrive at the result [7, 37] 



e T cot — = -Lm cot — 

27 4 27 



Z7 a 

where we have used (2.4) and (2.33). 

We now consider the boundary energy. From (B.5), we have 



(2.46) 



e-'^LB(^'{k) + e'^LBW'{k) 



i-zk) 

7(W++1) 

+e 2 



^7(«++l) 



7(ff_+l) 



-7r)fe 



+ e- 



7(H-+1). 



(2.47) 
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Substituting this contribution from (2.44) into (2.38), we obtain the boundary energy 



fOO 

— I dk e 

la 



9 



-iAk / sinh 



■yk 



1^ cosh ^ sinh 



Trfc 
2 



cosh ( (7(//+ + 1) - 7r)| ) + (//+ ^ //_) 



+ 



sinh ^ cosh((7 - f )A;) 1 



cosh ^ sinh ^ 



cosh 



7A; 



(2.48) 



which can be evaluated using the same contour integral as before, 



m 



1 + cot — + 
47 



^2 sin f(if+-^) sin f(i/_-^ 



+ 



sm 



27 



sm 



27 



(2.49) 



This matches with the result in [9] . 

The Casimir energy is given by (see (2.38), (2.44)) 



g dk 
En^-- ^ e 



a 



-iAk 



1 



27r cosh 



'yk 



LA'{k) + CA'{k) 



(2.50) 



Passing to coordinate space and taking the continuum limit, we obtain 



E, 



c 



2£ 

07 

m 
27' 



/>( 



1 



cosh - (A — x + ie) 



In A{x — ie) 



dxe^^^^''hnA{x-ie) = d9 e^~'nnA{9 - ie) , {2.51) 



27r 



where the renormalized rapidity 9 is given by (2.34). Using A{x) = A(— x) and ^mz 
—^mz, the last expression can be rewritten as 



m 



Ec = --^m / d9 \e^-'' In A(^ - ie) - e-^+'' In A(-^ + ie)] 



m 
2^ 



/oo 
d9 smk{9 -ie)hik{9 -ie) 
-00 



(2.52) 



3 Spin-1 XXZ / SSG with Dirichlet boundary conditions 

We turn now to our main interest, the spin-1 XXZ/supersymmetric sine-Gordon model with 
Dirichlet boundary conditions. 
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3.1 T — Q equations 



For the spin-1 chain, there are two relevant commuting transfer matrices: Ti{x) with a spin- 
I (two-dimensional) auxiliary space, and T2{x) with a spin-1 (three-dimensional) auxiliary 
space. The corresponding eigenvalues (which we denote by the same notation) obey T — Q 
equations found in [35]: Ti{x) can be written as ^ 

T,{x) = h{x)+l2{x), (3.1) 

Ux) = smh(2x + zj)B^-\xU(x + , 

Q{x) 

hix) — sinh{2x — i^)B^~^\x)(l){x — i^) — . , — , 

Q[x) 

and T2{x) can be written as 
T^ix) = X,{x) + X^ix) + Xs{x) , (3.2) 

X,{x) = sinh(2a; - 2t^)B^+\x - '^)B^'-\x + ^)0(a; - ^)0(a; - t) ' 

X2{x) = sinh(2x)S(-)(x-^)5W(x + 5^)0(x-5^)0(x+^)^^" + '^^^^""'^^ 



A3(x) = sinh(2x + 2i7)sH(x-^)sH(x+5^)0(x + ^)0(x+^^^^^-^ 



where 



0(a;) = sinh^(x- A)sinh^(x + A), 
B'^^\x) — sinh(x ± i?7+) sinh(x ± i77_) , (3.3) 

M 

Q{^) = W sinh(a; — Vk) sinh(a; -|- Vk) ■ 

k=l 

We denote the bulk and boundary parameters by 7 and r]±, respectively; A is the inho- 
mogeneity parameter which provides a mass scale; N is the number of spins; and the M 
zeros Vk of Q{x) are the Bethe roots. As we shall see below (4.4), it suffices to restrict 7 
to the domain < 7 < |. The domains (0 , |) and (f , f ) correspond to "repulsive" and 
"attractive" regimes of the SSG model, respectively. Note that T2{—x) = —T2{x). 

The fusion relation 

Ti(x - ^) T,{x + ^) = fix) + To{x) T2{x) , (3.4) 



^In [35], the transfer matrices were defined with some multiplicative factors which we omit here. 
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where 

To{x) = sinh(2a;) , 

fix) = l,{a:-'-l)h{x+'-l) (3.5) 

= sinh(2x - 2ij) sinh(2x + 2i-f)B^+\x - ^)B^-^ {x + ^)(P{x - ^)(p{x + ^) , 
can be readily verified using the identities 

A2(x) = + (3.6) 

For the homogeneous case A = 0, the local Hamiltonian (1.2) is obtained from the first 
derivative of the transfer matrix T^ix) [35]. For the inhomogeneous case A 7^ 0, as in the 
spin- 1/2 case (2.3), we define the energy by 



(3.7) 



where a is the lattice spacing, and g is given by (2.4). 

We consider the ground state. For simplicity, we restrict the boundary parameters ?7± to 
the interval 

|<77±<7r. (3.8) 

We argue in Appendix A. 2 that the boundary parameters should be further restricted to the 
range 

y + 7<77+ + 77-<Y + 7, (3.9) 

in order for the ground state to have M — N Bethe roots Vk which form approximate 
two-strings; that is, pairs Xk ± iyk with real centers Xk and imaginary parts i/k satisfying 
0<?/fc — ■^<Sl,as shown on the left-hand side of Figure 3. 

Since Q{x) can have zeros near the real axis (namely, when 7 is close to zero), it is 
convenient to work instead with the shifted quantity q{x) defined by 

q{x)^Q(^x+ '^^ , (3.10) 
14 



QM q(x) 
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- i3t/2 








-iir/2 



Figure 3: Schematic depiction of zeros of Q(x) and q{x), respectively 

which is ANZ near the real axis. Also, one can check numerically that Ti{x) and T2{x) do 
not have zeros near the real axis except for a simple zero at the origin. To remove this root, 
we define 

f.W = Ilfel, f,W = ZM, (3.11) 

j2[x) jdyx) 

where again is any function whose only real root is a simple zero at the origin, in 

particular /i(0) = 0, /x'(0) ^ 0, so that Ti(x) and ^2(0:) are ANZ when x is near the real 
axis. 



3.2 Auxiliary functions 

It is convenient to define the auxiliary functions [25] 

Ai(x) + A2(x) - A3(x) + A2(x) 

= VT~\ ' Kx) = = K~x)- 3.12 

Asia;) Ai(x) 

Since h{x) is the complex conjugate of h{x) for real x, we shall generally refrain from writing 
equations for 6(x), as they can be readily obtained by conjugation of corresponding equations 
for h{x). With the help of (3.6), we obtain 

Ti{x-'i) (l>{x-'i) B^+){x + 'i) g(x + f^) 

sinh(2x + 2*7) 0(x + |)0(x + f^) Bi-){x - ^)5(-)(x + |) Q{x - f^) ' 

In terms of the quantities (3.10), (3.11), we rewrite this relation in the compact form 



h{x) = C,{x) Ux - ^) ^|^^|^^j , (3.14) 
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with 

' sinh(2a; + 2i-i)<p{x + f + f^) B(-) {x - (a; + |) ' ^ ' ' 

Defining 

= l + = 1 + = , (3.16) 

we also obtain 

T2{x) = Asia;) B{x) = \i{x) B{x) . (3.17) 
In terms of the quantities (3.10), (3.11), we can reexpress (3.17) as 

Ux) = Ux)'j^^-^^-^B{x), (3.18) 
'-y^^r-^Bix) , (3.19) 



with 



We note that i?(x) has zeros just above the real axis, and B{x) has zeros just below the 
real axis. Indeed, due to the factor Q{x + y) in the definition of X3{x) (3.2), and the fact 
that the imaginary parts of the Bethe roots have magnitude > ^, we see that Xs^x) has 
poles just above the real axis (and also just below the line '^mx = —7). It follows from 
(3.17) that B{x) must have corresponding zeros to cancel these poles, since the product 
A3(a;)-B(a;) = T2{x) is analytic. Similarly, since Ai(a;) has poles just below the real axis (and 
also just above the line '^mx = 7), B{x) has corresponding zeros to cancel these poles. 

Finally, we define the auxiliary functions y{x) and Y{x), 

^ T,{x)Ux) Y{x)^l + y{x), (3.21) 

in terms of which the fusion relation (3.4) becomes 

^/ ^7n^/ il s f(x)Y(x) ,„ „, 

T,{x - -1) T,{x + ^) = , '\' . (3.22) 

2 2 - ^) n{x + ^) 
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3.3 Derivation of NLIE 



Since lnT2(x) is analytic near the real axis, Cauchy's theorem gives 



0= 6 dx [\nT2{x)]"e"''' = / [\nT2{x)]"e"'''dx + / [\nT2{x)]"e"'''dx , 



(3.23) 



where we choose the contour C as in Figure 2, with e such that max [yk] — ^ < e. These 
integrals can be written using (3.18) and (3.19) as 



/ [hif2{x)Ye'^''dx = f dx {lnt_(a;)}"e^^^+ / da; i In 
Jci Jci Jci I 



+ / dx {lnB{x)}" e''^ 

'Ci 



IC2 



'C2 



/„'■ 

















g(x + f -f) 



ikx 



(3.24) 



g(x+3f -f) 



Akx 



+ 



/ dx {lnS(a;)}"e^'=^, 

JC2 



(3.25) 



so that the quantities q{x + y — y) and q{x — ^ + y) are integrated along Ci and C2, 
respectively, and not the other way around. ^ 

In terms of Fourier transforms defined by (2.15), we can rewrite (3.24) and (3.25) as ^ 

e(i-f)^-e(^-^)H 



uFll{k) = LiJ(A;) + 



(^-f)^_p(f-i)^ 



Lq"{k) + CB"{k) , 
L^'{k) + LB"{k) , 



(3.26) 
(3.27) 



respectively. In obtaining these results, we use the periodicity of q{x) as in (2.17) to make 
the imaginary part of the argument negative. 

Adding (3.26) and (3.27), and remembering the fact (3.23) 



£f^'(fc) + LT2"(^)=0, 
we obtain an expression for Lq"{k), 



(3.28) 



A.2 2 — q\2 2)'^ — q\2 2)'^-\-q\2 2)" 



Lq"{k) = LB"{k) + CB"{k) + Dq{k) , (3.29) 



^Because of our choice of e, the quantities q{x+^ — ^) and B{x) are nonzero along Ci, and q{x— ^ + ^) 
and B{x) are nonzero along C2. 

^Note that all the integrals of q have been expressed in terms of Lq"(k). This would not have been 
possible if we had interchanged Ci and C2 in (3.24) and (3.25). 
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where we define 



Dg{k) = Lt'[{k) + Ct'L{k) 



(3.30) 



Prom the expression (3.14) for b{x) and the fact that Ti{x) is ANZ near the real axis, we 
obtain 



Lb"{k) = e--LTi'{k) + Lq"{k) 



.(^-f)^_p(f-¥)^' 



+Lci:{k) 



(3.31) 



Prom the fusion relation (3.22), we obtain 

(e'r + e-^ ) L^l'{k) = LV"{k) + Lp'{k) - + e'"^^ Li?'{k) + e^5{k) , (3.32) 

where 5{k) is again given by (B.3). Substituting Eqs.(3.29) and (3.32) into Eq.(3.31), we 
obtain 



Lh"{k) = 



e(¥-f)fc_e(f-¥)^ 



l\2 2 — (>\2 2)'^ — f>\2 2)'^-\-(A.2 2 



+ 



Tfc 

e 2 



7fc Tfc 

e 2 + e 2 



LS"(A;) +£S"(A;) + D5(A;)J 



LY"{k) + L/"(/c) + Lq'(A;) - e-'^L^"{k) + 



7fc 2h. 

e 2 + e 2 



where 



G'(^) 
G2{k) 



-G{k) LB"{k)+CB"{k) +G2{k)LY"{k) + C{k), 

sinh ((7r-37)|) 
2 cosh fsinh ((71-27)1) ' 



e 2 



e 2 + e 2 



(3.33) 

(3.34) 
(3.35) 



C{k) = 

We find that C{k) is given by (B.ll 



G{k) D,{k) + G2(A;) L/"(^) + LC'^{k) - e-'^Li,"{k) + J^^\^ .(3.36) 



e 2 + e 2 



2-Kk< N 



2 cosh f 



+ 



[sinh ((r/+ - f )A;) + sinh ((ry_ - ^)k)] 
2coshf sinh((f -7)A;) 



+ 



cosh ^ sinh ((37-7r)|) 
cosh f sinh ((27-7r)|) 



(3.37) 



We now turn to the third and final NLIE equation. Prom the definitions of y{x) (3.21) 
and T2{x) (3.11), we obtain 



Ly"{k) = Lf^'ik) + LT^'ik) + Lfi"{k) - Lf"{k) . 



(3.38) 



18 



To evaluate LT2{k), we combine Eqs. (3.26) and (3.27) to cancel Lq"{k), namely, 

which together with (3.28) gives 

(|e^ + ) Lf^'ik) = e-'-^LB"{k) - e^CB"{k) + Dr^k) , (3.39) 

where we define 

D^ik) = e-'^Ltlik) - e^C?{k) . (3.40) 
Substituting the result (3.39) for LTll^k) into (3.38), we obtain 

L^'{k) = -G2{-k) CB"{k) + G2{k) LB"{k) + Cy{k) , (3.41) 



where we define 



We find (B.13) 



Cy{k) = -^^^ + LTHik) + L^,"{k) - Lf"{k) . (3.42) 

e 2 + e 2 



Cy{k) = AT:kG2{-k) . (3.43) 



In summary, the NLIEs of the lattice SSG model with Dirichlet boundary conditions in 
Fourier space are 

Ly'{k) = -G{k) \L§"{k) + CB"{k) \ + G2{k) LY"{k) + C{k) , (3.44) 

li'^k) = G2{k)LB"{k)-G2{-k)£B"{k) + Cy{k), (3.45) 

where C{k) and Cy{k) are given by Eqs. (3.37) and (3.43), respectively. Passing to coordinate 
space, integrating twice, and taking the continuum limit, we obtain 

/OO /"OO 
d9' G{9 -9' - is) In B(^' + ie)- / d9' G{9 -9' + ie) In B(^' - is) 
■OO J — oc 

/OO 
d9' G2{9 -9' + ie) In Y(^' - ie) + i2mL sinh 9 + i Pbdry{9) - in , 
-OO 

d9' G{9 -9' - ie) In B(^' + ie) + / d9' G{9 -9' + ie) In B{9' - ie) 

-OO J —03 

/OO 
d9' G2{9' -9 + ie) \nY{9' + ie) - i2mLsmh9- i Pbdry{9) + in , 
-OO 

/OO /'OO 
d9' G2{9 -9' + ie) In B{9' - ie) + / d9' G2{9' -9 + ie) In B{9' + ie) 
-OO J — OO 

+ ^Pj;(^). (3.46) 
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As in the spin- 1/2 case, the continuum hmit consists of taking A — > oo, N ^ oo and lattice 
spacing a — > 0, such that the interval length L = x+ — and the soliton mass m are 
given by (2.33). The renormalized rapidity 9 is again given by (2.34), and we have defined 
b(^) = b{^), etc. as in (2.36). Hence, the kernel G(^) is given by 

G{e) = / dk e-^*^^^/'^ d{k) , (3.47) 
where G{k) is given by (3.34); and G2{0) is defined similarly in terms of G2{k) (3.35), 

where 6 has a slightly positive imaginary part. Finally, Phdryid) is given by 

Pbdry{0) = I ^ dx'R{x') = - dx'R{x') = ^ I dO' I dk e-'^^^'l^ R{k) , (3.49) 

TT 

where R{k) is given by 



J [sinh ((77+ - l)k) + sinh ((ry- - f )A;)] cosh ^ sinh ((37 - 7r)f) 
\ 2coshf sinh ((77-27) I) ^ coshf sinh((27-7r)|) 



and Py{0) is given by 

i-e a 

Py{e) = 47r d^' G2(-^') = -2«lntanh- -27r, (3.51) 

J-oo 2 

where 9 has a slightly negative imaginary part. The integration constants are explained in 
Section C.2. 

These NLIE equations are similar to those for the periodic chain [23, 25, 28], with addi- 
tional boundary terms involving Pbdry{9) or Py{0). These terms constitute one of our main 
results. As we shall see, these boundary terms make essential contributions to the boundary 
S matrix (IR limit) and to the effective central charge (UV limit). 



3.4 Vacuum and Casimir energies 

The energy computation is similar to the one for the spin-1/2 case in Section 2.3. We see 
from (3.21) that T2{x) can be expressed in terms of y{x), 

T2{x)-^^r^. (3.52) 
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It follows from the energy definition (3.7) together with the fact (2.39) that 
Recalling the results for Lf"{k) (B.9) and Ly"{k) (3.45), we obtain 



a I 27r 



g-ifc(A+^) _ g-jfc(A-rt) 



(3.53) 



E 



a. 27r 



+ 2 sinh ( — 



e-"-T LBi+y {k) + e'TLBi-Yik) + 4mG2{-k) 

1 



[—ik) ^ ' 



cosh2i 



LB'{k)+ LB'ik) 



(3.54) 



The first term gives the bulk vacuum energy, which can be written explicitly as 

2Ng 



la 



-oo 



sinh ^ 



(3.55) 



The second term gives (with the help of (B.IO)) the boundary vacuum energy 

1 



Eu 



ia 



dke ''^''sinh 



(t) I ^ (^^^ - i - + ^ ^-l 



+ 



7fc 

e 2 



coshf 



+ 



;(f-7)fc + e(7-|)fc _ 
2 sinh ^ 



(3.56) 



and the third term gives the Casimir energy. 



9 



En^-- / — e 



dk 



-iAk 



.ooSTr"- coshf 



LB'{k) + CB'ik)] . 



(3.57) 



In order to take the continuum limit, we adopt (as in the spin-1/2 case) the renormal- 
ization procedure of keeping only the (finite) terms that can be expressed in terms of the 
physical mass m (2.33). We implement this procedure by closing the integral contours in the 
lower half plane, and selecting only the contribution from the residue at A; = — Since the 
integrand in (3.55) is analytic at k — — — , the bulk energy vanishes. 
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Eb = 0, 



(3.58) 



in agreement with known results (see, e.g., the second reference in [23]). We obtain from 
(3.56) that the boundary vacuum energy does not vanish, however, and is given by 



Eh = m . 



(3.59) 
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That is, each boundary contributes the vacuum energy m/2. Note that this result is in- 
dependent of the boundary parameters. We shall present further support for this result in 
Section 6. Finally, we find that the Casimir energy (3.57) is given by 

Ec = —^m / de sinh(^ - ie) In B(^ - is) . (3.60) 



4 Infrared limit 



In this Section we analyze the IR limit mL — > oo. In particular, we compute the SSG 
soliton boundary S matrix (4.27), and show that it coincides with one proposed by Bajnok 
et al. [32]. In making this identification, we determine the "lattice <-> IR" relation for the 
boundary SSG parameters (4.25). 

Before starting this computation, we note the relations among the three hulk SSG pa- 
rameters (see Figure 1). Let A denote the IR bulk SSG parameter, and recall that j3 and 7 
are our UV and lattice bulk SSG parameters, respectively. (See Eqs. (1.1), (1-3).) The SSG 
soliton bulk S matrix [22] has the product form SG(9\\) RS0S{9), where SG{6; \) is 
the SG soliton bulk S matrix [1], and RS0S{9) is the quantum-group restricted SG model 
bulk S matrix [2]. The bulk "UV ^ IR" relation is given by [22], 



The bulk "lattice <-> IR" relation is given by 

1 



E - 2 

7 



A. (4.2) 



This relation can be readily inferred from a comparison of the kernel G{k) (3.34) with the 
integral representation of the SG soliton-soliton scattering amplitude SG++{9 ; A), 

i de ^ ' - y_ ' 2 cosh f sinh f ■ ^^-^^ 

The relations (4.1), (4.2) imply the bulk "UV ^ lattice" relation 

= 4 (tt - 27) . (4.4) 

The condition < < 47r therefore implies < 7 < |. The free Fermion point A = 1 
corresponds to 7 = |. In the "repulsive" domain 7 G (0 , |), the particle spectrum consists 
only of supersymmetric multiplets of solitons and antisolitons. In the "attractive" domain 
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7 e (f ) f )) the spectrum also includes bound states, namely, supersymmetric multiplets of 
breathers of mass [22] 

mn = 2msin(^— j, n = 1 , 2 , . . . , [AJ . (4.5) 

We turn now to the boundary theory. Let ^± denote the IR boundary SSG parameters, 
and recall that ip± and r]± are our UV and lattice boundary SSG parameters, respectively. 
(See Eqs. (1.6), (1.5).) 

We define the boundary S matrices (refiection factors) R{9h ; A , for a soliton with 
mass m and rapidity 6h by the Yang equation for a single soliton on an interval of length 
L > 1/m, 



e 



i2mLsin^e,j^^g^ ; A , C-) R{9h X , ^+) ^ 1 ■ (4.6) 



We shall compute these S matrices by deriving a similar relation from the IR limit mL — > oo 
of the NLIE for a state of one hole with rapidity 9h, 

In b(^) = i2mL sinh 9 + iPbdryi9) + ixi9 - 9h) + ix{9 + 9h) 

/oo 

d9' G2{9 -9' + is) In Y(^' - is) - iir , (4.7) 
-oo 

Iny(^) = tPy{e) + zgy{9-9h)+tgy{9 + 9h), (4.8) 
where x{9) and gy{9) are the hole source terms [25, 28] 

X{9) = 27r d9' G{9') , gy{9) = -i In tanh ^ + | , (4.9) 

and where, in the latter equation, 9 has a slightly negative imaginary part. Indeed, since 
lnh{9h) is in times an odd integer (see, e.g., [25, 28]), evaluating Eq. (4.7) at 9h and 
exponentiating both sides gives 



e 



i2mLsmheh ^iPbdry{6h)+ix('^6h)+)C(eh) _^ 



where }C{9) is the convolution term in (4.7), 

/oo 
d9' G2{9- 9' + is) InY {9' -ie). (4.11) 
-oo 

Comparing (4.10) with the Yang equation (4.6), we conclude that the product of boundary 
S matrices is given by 

R{9h ; A , C-) Ri9h ; A , C+) = e'^'^^y^^^^+'''^'"^^^^+'^^'^'^^ . (4.12) 
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We evaluate first the factor e'^^^''^ which, as we shall see, is the RSOS factor. To this end, 
we observe from (3.51), (4.8) and (4.9) that y(^) is given by 



o9 , 1 1 



y{e) = - tanh^ - tanh -{6- 9^) tanh -(^ + Oh) . 



It follows that 

InY(^) =ln(l+y(e)) = ln 



cosh Y cosh 6 



cosh^ f cosh Ue- Oh) cosh + Oh) 



(4.13) 



(4.14) 



The convolution term (4.11) is therefore given by 



]C{e) = In cosh ^ + Qi{e) - 2Q2{e) - Q2{0 - dh) -Q2{e + en) , 



(4.15) 



where [28] 

Qi{0) 
Q2{e) 

with 



i 1 
d9' G2{9 -9' + ie) In cosh(e' - ie) = -- tan"^ sinh ^ + - In cosh 9 , 

d9' G2i9-9' + ie)lacoshh9' -ie)^~X2{9) + \lncosh^, (4.16) 



dk e 



-ike 



4cosh=^^ ' 



X2(0) = l 



It follows that 



^Oh) = - [X2{29h)+2x2{9h)-t&^-'smh9n] 



(4.17) 



(4.18) 



Differentiating with respect to 9h, and then making use of (4.17) and the Fourier transform 
result 

cosh 9 



we obtain 

d 



2x'2{29h) + 2x'2{9h) 



1 



cosh 9h 



2 cosh ^ ' 



(4.19) 



-2ik0h 



dk 



cosh^ ^ cosh^ Ttk 



Upon integrating, we conclude that 



min\"h) i 



(4.20) 



(4.21) 



^We use ~ to denote equality up to crossing factors of the form e'^°"'** ^. Such factors have also not been 
obtained in the bulk case [28]. 
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where Pmin{&) is a reflection factor of the boundary tricritical Ising model [31], whose integral 
representation is given by [39] 



,11 " dt sm(2t9/7c) 



(4.22) 



We stress that the boundary term Py{9) in (4.8) is essential for obtaining this result. Since 
Py{0) is independent of the boundary parameters, so is the result (4.21). 

We turn now to the remaining factor e^^^''-''^^^^^'^^^^'^^^^ in (4.12) which, as we shall see, is 
the sine-Gordon factor. Differentiating with respect to 9h, and recaUing the Fourier transform 
results (3.34), (3.50), we obtain 

^ [PMryien) + xi^Oh)] = P'MryiOH) + 2x'(2^,) 



dke' 



sinh - f )^ + siiili - f )- 



2 cosh f sinh ((f -2)f ) 



+2 



sinh ^ sinh ((^ -3)^)) 



sinh nk sinh 



((f-2)f 



(4.23) 



Let us compare this result with the soliton reflection amplitude P+{9 of the boundary 
SG model with Dirichlet boundary conditions [3] , which has the integral representation [38] 



Id r°° 
I de 



sinh ((1 + f|)^) sinh ^ sinh ((i - 1)^)) 



2 cosh f sinh g 



+ 



sinh Tik sinh ^ 



, (4.24) 



where A and ^ are the bulk and boundary IR parameters, respectively. Recalling the bulk 
"lattice <-> IR" relation (4.2), and assuming the boundary "lattice <-> IR" relation 



{-?) 



(4.25) 



we conclude that 



(4.26) 



Combining the results (4.12), (4.21) and (4.26), we conclude that the NLIE generates the 
following SSG sohton boundary S matrices 



(4.27) 
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This is the boundary S matrix which was proposed by Bajnok et al. [32] for the the Dirichlet 
BSSG"*" model. This is another of our main results. Similarly to the bulk case, the SSG 
boundary S matrix is a product of SG and RSOS boundary S matrices. 

For general values of boundary parameters, the BSSG+ model with one boundary has 
the conserved supercharge [30, 31, 32] 

g+ = g + Q + yr , (4.28) 

where V = (—1)'^ is the Fermionic parity operator, and 7' is an undetermined parameter. ^ 
Bajnok et al. also propose the relation 

Q\^2(^H + mZ^ , (4.29) 

where H is the Hamiltonian, and Z is the topological charge. Since the ground state has 
g+ = 7' and Z = 0, it has energy H = 7'^/2. Our result (3.59) that each boundary 
contributes vacuum energy m/2 implies that 

7' = ±v^, (4.30) 

at least for the Dirichlet case. That is, we have succeeded to fix the undetermined parameter 
in the scattering theory proposed in [32]. 



5 Ultraviolet limit 

In this Section we first analytically compute the Casimir energy Ec{L) (3.60) in the UV 
limit mL — > 0. The result (5.17) is proportional to the effective central charge 

^c(0) = -^c,//(0), (5.1) 

where Ce// (0) = c — 24Ao, c is the central charge, and Aq is the Lq eigenvalue of the ground 
state. We then compare this result for Ce//(0) to the value for the conformal limit of the 
SSG model with Dirichlet boundary conditions. In this way, we obtain a boundary "UV ^ 
lattice" relation (5.22). When combined with the boundary "lattice IR" relation from 
the previous section (4.25), we obtain the boundary "UV IR" relation (5.23). 

^This parameter is called 7 in [32]; however, here we add a prime in order to distinguish it from our bulk 
lattice parameter. 
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5.1 NLIE computation 



We now proceed to analytically evaluate the Casimir energy in the UV limit mL — > 0. As is 
well known, only large values of |^| contribute in this limit. Let us first consider 9^1, and 
define the finite rapidity 9 by 



9^9- ln{mL) . 

The corresponding contribution to the Casimir energy (3.60) is given by 



LE. 



c 



IGtt 



d92iJ 



lnB+(^') -lnB+(^) 



(5.2) 



(5.3) 



where the auxihary functions B+(^) = B(^ — In(mL)), etc. satisfy the NLIE equations 

/OO POO 
d9' G{9 -9' - ie) In B+(^' + ie)- / d9' G{9 -9' + ie) In B+(^' - ie) 
-OO J —OO 

/OO 
d9' G2{9 -9' + ie) In Y+(^' - ie) + ie^ + iPbdry{oo) - in , 
-OO 

/OO roc 
d9' G{9-9' - ie) In B+(^' + ie) + / d9' G{9 - 9' + ie) In B+(^' - ie) 
-OO J —OO 

/OO 
d9' G2{9' -9 + ie) lnY+(^' + ie) - ie^ - iPbdryioo) + in, (5.4) 
-OO 

/OO poo 
d9' G2{9 -9' + ie) In B(^' - ie) + / d9' G2{9' -9 + ie) In B(^' + ie) , 
-OO J —OO 

which are independent of mL. The P^^ term (3.51) docs not appear in the third equation 
due to the fact Py{oo) = mod 2n. We note here for later reference that 



1 vr 

Pbdry{oo) = -R{0) = — {r]+ + r?_ - 37) . 

z n — 



(5.5) 



We use the so-called dilogarithm trick [6, 24]. We first rewrite the NLIE equations (5.4) 
in matrix form, 



u — V + K * w , 



(5.6) 



where 





( In b+ ^ 




/ 


ie^ + iPbdry{oo) - in 




/ lnB+ 






u — 


lnb+ 


I- 




—ie^ — iPbdryi^) + 


, w = 


lnB+ 




(5.7) 




V iny+ ) 




\ 


) 






) 
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the kernel K is symmetric (i.e., Kij{9 , 9') — Kji{9' ,9)), and the star * denotes convolution. 
We see from (5.6) that ^ 



T / /T \ 

[W V — W Vj , 



d9 {w^u' - w'^u) = d9{i 

1 J —oo 

since the symmetry of the kernel implies that 

/oo 
S {w'^K' *w- w'^K * w) = . 
-oo 

It follows from (5.8) that 

I d9i (In b+)' In B+ - In b+(ln B+)' + (In b+)' In B+ - In b+(ln B+)' 

J —oo 

+(lny+)'lnY+ - lny+(lnY+)'| = f d9 2ie^ (in B+ - In B+) 

J —oo 

-i (e^ + Pbdryioo) - tt) (in B+ - In B+) 



(5.8) 



(5.9) 



e=oo 
e=-oo 



(5.10) 



Since B+ = 1 + b+, etc., the LHS of (5.10) can be expressed in terms of the dilogarithm 
function L+{x), defined by 



-ln(l + y) - —^\ny 

y i + y . 



(5.11) 



The integral on the RHS of (5.10) is essentially the sought-after quantity LE^ (5.3). We 
conclude that 



LE. 



c 



= ^{L+(b+(«D))-L+(b+(-oo)) + L+(b+(«D))-L+(b+(-oo)) (5.12) 

-OO J 



+ L+ (y+(oo)) - L+ (y+(-oo)) + 2 + ^bdry{oo) - ttJ (in B+ - In B+) 



The plateau values of the auxihary functions can be obtained from the NLIE equations 
(5.4). For ^ — > 00, we readily obtain 



b+(oo) = b+(oo) — , y+(oo) — 1 ■ 



(5.13) 



To determine the plateau values for ^ — — cxd is less trivial, as the corresponding plateau 
equations are nonlinear. We make the Ansatz 

b+(-oo) = e^'^ (1 + e"^) , b+(-oo) = e'"^ (l + e"^'^) , y+(-oo) = 1 + e^'^ + e'"^ ,(5.14) 



^Here the prime denotes differentiation with respect to 6. 
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where u is still to be determined. We indeed find a solution with 

7-7r) , 



P6drj/(oo) - TT 

OJ = — i--^^ — - — - = 2 (77+ + rj. 



I-G(O) 



(5.15) 



where we have used the result (5.5). For the plateau values (5.14), the following sum rule 
holds [25] 

^2 



L+ (b+(-oo)) + L+ (b+(-oo)) + L+ (y+(-oo)) - L+(l) 



TT 



m < 



27r 



(5.16) 



The above bound is satisfied when the boundary parameters are in the domain (3.9). Noting 
also that -L+(0) = 0, we conclude that the Casimir energy is given by 



TT 



LEc{Q)^2LE^^-- 



12 



2 7r(7r - 27) 



(77+ + 77_ - 7 - T^y 



(5.17) 



In obtaining the first equality, we have used the fact that the contribution Eq from ^ <C — 1 
is the same as E'^, and that -Bc(O) = E'^ + E^. 



5.2 CFT analysis and UV-IR relation 

Our result for the Casimir energy (5.17) together with the relation (5.1) evidently imply that 
the effective central charge has the value 

Ce//(0) = c - 24Ao = \- ^^^l^^^ {v+ + V--^- ■ (5.18) 

We remark that for boundary parameter values 77+ — rj^ — |(7r + 7), the boundary terms 
in the Hamiltonian (1.5) which are proportional to 5"^ vanish; and also Aq vanishes, so that 
Ce//(0) — c — 3/2. A similar phenomenon was observed for the spin-1/2 case in [9]. 

In the UV limit, the boundary SSG model with Dirichlet boundary conditions (1.1), 
(1.6) evidently reduces to a system of one free Boson and one free Majorana Fermion, each 
with Dirichlet boundary conditions. For the former model, the central charge and lowest 
dimension are given by [40] 

Cb = 1 , Ab = i- ((^_ - <f+f , (5.19) 

while for the latter model 

CF = ^, Ap = 0. (5.20) 

Since both the left and right boundaries have the same (Dirichlet) boundary condition, the operator 
content includes the identity operator, which has dimension zero. 
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It follows that the SSG model with Dirichlet boundary conditions should have 

c^cb + cf^^, Ao = Ab + Aj. = ((^_ - ip+f . (5.21) 

Comparing this CFT result with the NLIE result (5.18) and recalling the bulk "UV lattice" 
relation (4.4), we obtain the boundary "UV lattice" relation 

V±^l(7r + ^)±^^^. (5.22) 

Combining this result with the boundary "lattice ^ IR" relation (4.25), we finally arrive at 
the SSG boundary "UV ^ IR" relation 

i± = ±j^±- (5.23) 

This is another of our main results. 

The relation (5.23) is similar to the one found by Ghoshal and Zamolodchikov [3] for the 
SG model (namely, ^ = ^^Po)^ ^-nd it can be understood in a similar way. Indeed, for the 
SSG model, it is also plausible to assume a linear relation between these parameters. 



a + b(fi±. (5.24) 



When (/?-(- = 0, the model has the symmetry (/? i— > — (/?; and, since the RSOS factor of 
the boundary S matrix is proportional to the identity for BSSG+ [32], the soliton and 
antisoliton reflection amplitudes should be equal, which corresponds to ^± = 0. Thus, a = 0. 
Furthermore, as in the SG model, there arc boundary bound states corresponding to poles 
of the boundary S matrix at 6' = iun, where 

.„ = |-^^^, » = 0,1,.... (6.25) 
Moreover, these states satisfy [32] 

.2 . fl" 



= 2 y-^ + m cos Un + mj , (5.26) 

which is consistent with (4.29), since these states have Z = 1. When ip± have the half-period 
values ip± = ±7r//9 the n = bound state should have the same (5+ eigenvalue as the 
ground state. From (5.26), we see that this condition corresponds to z/q = ^! which in turn 
implies {± = 7r(A + ^) = 27r^//3^, as follows from (5.25) and the bulk UV-IR relation (4.1). 
It follows that b — ±27r//3, and so we recover the boundary UV-IR relation (5.23). 



^^The Lagrangian (1.1) has the periodicity if ip + 27r//3. 
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6 Intermediate volume and BCPT 

For intermediate values of volume / = mL, the Casimir energy cannot be computed analyt- 
ically. Nevertheless, it is possible to solve the NLIEs (3.46) numerically by iteration, and 
evaluate the Casimir energy through (3.60). Two sample plots of Ce//(/) = —2ALEc{L) /tx vs. 
In / are shown in Figure 4. The numerical result for Ceff{l) in the UV region / coincides 
with the analytical result (5.18). Also, as expected, Ceff{l) decreases monotonically to as 
I varies from the UV region to the IR region / — oo. 




-3 -2.5 -2 -1.5 -1 -0.5 0.5 



Figure 4: Ce// vs. InZ, for 7 = 27r/7, 77+ = !]_ = (vr + 7)/2 (solid 
line) and 77+ = ?7_ = 7r/2 + 37/4 (dotted line) 

For small values of /, these numerical results can be compared with those from bound- 
ary conformal perturbation theory (BCPT). We shall follow closely the presentation in the 
Appendix of the second reference in [11]. We regard the SSG model (1.1) as a perturbed 
boundary CFT, 

L = LbCFT + Lpertb , Lpertb = ^ / dx^{x,t) , (6.1) 

^ Jo 

where Lbcft is the Lagrangian for a free scalar field ip and a free Majorana Fermion field \1/ 
obeying Dirichlet boundary conditions (1.6). We restrict our attention here to the particular 
case V9_ = = (y9o; Order to avoid introducing boundary-changing operators. The bulk 
perturbing operator $(x, t) = cos{l3ip)^'^ is a primary field with weights (A, A), where 

1 /?2 1 

A = - + ^, -<A<1. (6.2) 

2 Stt' 2 ^ ' 
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The corresponding Hamiltonian is 

H{L) = Hbcft ■ ^^-^^ 

We map the infinitely long strip of width L to the upper half plane by setting z = e*^(^+*\ 
where t = —iy is the Euclidean time. Taking the Hamiltonian at t = and changing the 
integration variable to 6* = jx, we have 

with c — 3/2. First-order perturbation theory implies that the ground-state energy is given 

by 

Eo{L) = ^ (Ao- ^) + ^ {£l"^~" J\9{0me^',e-^<^m + O{ml), (6.5) 

where |0) is the ground state of the unperturbed theory. Using the fact that the bulk one- 
point function has the form 

mie'',e~^'m^j^^^, (6.6) 
and performing the 9 integration, we obtain 

where Ce//(0) — c — 24Ao, as before. The energy Eq{L) is the sum of bulk, boundary and 
Casimir energies 

Eo{L) = EbL + E, + Ec{L) , Ec{L) = -^c,ff{l) . (6.8) 

Recalling the bulk and boundary energy results (3.58), (3.59), we obtain 

24/ / TT \ 2A-2 Y{- - A^rfi) 

Ce//(/) = Ce//(0) + — - 3mo (-) r(l- A) """'^ + • ^^-^^ 

Following Cardy and Lewellen [42], the coefficient of the bulk one-point function can be 
expressed as a ratio of scalar products, 



Cbulk — 



{{0\D) {0,0\Bn{ipo)) 0=1 ((0|0) 
Here \Br){(po)) is the c—1 Dirichlet boundary state [40] 



. (6.10) 

Ising 



\BD{fo))^^fD J2 e^'''^''°e-S^=i^°-""-"|0,A;), (6.11) 



k=—oo 
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where \0,k) is annihilated by q;„ and q;„ for n > 0, and has weight /io,fe — ^g^; ^-nd is a 
normalization factor. Hence, e^^'^ (which have weight ^) are identified with |0,±1). Thus, 
I cos(/3(/7))) = I (|0, 1) + |0, —1)), and we conclude that 

{{cos{P<p)\BDi<Po)) 



{0,0\BdM) 



c=l 



= cos(/3(^o) ■ 



(6.12) 



Moreover, — —^e, where e is the Ising energy density operator. Recalling the expressions 
[43] for the boundary states corresponding to fixed boundary conditions in terms of Ishibashi 
states 



we obtain 



|0) 






1 




h 






((^|0) 






((0|0) 


Ising 











1 







cr 



(J 



1. 



(6.13) 



(6.14) 



We conclude from (6.10), (6.12) and (6.14) that the coefficient of the bulk one-point function 
is given by 



Cbuik = --cos(/3(/?o) ■ 

71 



(6.15) 



Moreover, the parameter mo is related to the SSG soliton mass m by the so-called mass-gap 
formula [41] 

8r(A) rvrm /2A-1\1'-'^ 

nil, where mi is the mass 



1-A J 

0, namely mo 



r(i - A) [ 8 

This relation has the correct classical limit (3 
of the first breather (4.5). 

Substituting (6.15) and (6.16) into (6.9), we obtain the result 

24 



Ce//(0 



Ce//(0) + ^ (I - ail 



2-2A 



where 



2A-1 

4(1^ 



2-2A r(i-A)r(i)r(A) 



cos(/?(^o) ■ 



(6.17) 



(6.18) 



r(i - A)2 

We have not attempted to compute the second-order correction a2- In Table 1, we compare 
the analytical BCPT result for cti (6.18) with results obtained by fitting numerical NLIE 
values for Ce//(Z) to the curve (6.17). The excellent agreement between the analytical and 
numerical values further supports the validity of our NLIE (3.46), as well as the boundary 
energy result (3.59) and the boundary "UV <-> lattice" relation (5.22). 



^^We evidently obtain the same result if we consider |i) instead of |0). 
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7r/7 




= C) 




= 0.1 


Y^O - 


= 0.2 


NLIE 


BCPT 


NLIE 


BCPT 


NLIE 


BCPT 


2.6 


1.05337 


1.05314 


1.03813 


1.0379 


0.992871 


0.992644 


2.8 


0.977854 


0.977706 


0.960355 


0.960207 


0.908481 


0.908335 


3 


0.898947 


0.898868 


0.880187 


0.880108 


0.824691 


0.82461 


3.5 


0.717426 


0.717416 


0.698197 


0.698184 


0.641543 


0.641519 


4 


0.572666 


0.572698 


0.554784 


0.5548 


0.502248 


0.502225 


4.5 


0.46177 


0.461973 


0.445794 


0.445941 


0.398955 


0.398957 


5 


0.37713 


0.377594 


0.363094 


0.363448 


0.321985 


0.322071 



Table 1: Comparison of NLIE (numerical) and BCPT (analytical) results for ai, 
for various values of bulk and boundary parameters. 

7 Discussion 

We have proposed a set of nonlinear integral equations (3.46) to describe the boundary 
supersymmetric sine-Gordon model BSSG"*" with Dirichlet boundary conditions on a finite 
interval (1.1), (1-6). In particular, we have found the boundary terms Pbdryi^) ^-nd Py{0) 
which encode boundary effects. We have computed the corresponding boundary S matrix 
(4.27), and found that it coincides with the one proposed by Bajnok et al. [32] for the 
Dirichlet BSSG"'" model, with conserved supercharge (4.28) depending on a parameter 7'. 
We have determined this parameter (4.30) by computing the boundary vacuum energy. We 
have also proposed a relation (5.23) between the (UV) parameters in the boundary conditions 
and the (IR) parameters in the boundary S matrix. Moreover, we have demonstrated that 
the NLlEs can be solved numerically for intermediate values of mL, and we have found 
agreement with our analytical result (5.18) for the effective central charge in the UV limit 
and with boundary conformal perturbation theory (6.17), (6.18). 

There are a number of related questions which remain to be addressed. While we have 
focused here primarily on the ground state, it should be interesting to study bulk (along the 
lines [28]) and also boundary excitations. It would also be interesting to formulate the TBA 
equations for this model, and compare the results with those from our NLIE. 

Ultimately, we would like to extend this (Dirichlet) analysis of BSSG"*" to the case of 
general integrable boundary conditions [30]. For this we would need the Bethe Ansatz 
solution of the open spin-1 XXZ quantum spin chain with general integrable boundary terms 
[34], which is currently under investigation. 

It is not evident which integrable spin-1 model, if any, corresponds to the Dirichlet BSSG~ 
model. Indeed, the boundary terms (1.5) already correspond to the most general diagonal 
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c-number solution of the boundary Yang-Baxter equation. Perhaps additional dynamical 
boundary degrees of freedom or impurities will be required. An interesting feature of BSSG~ 
is that, in contrast to BSSG"^, the RSOS part of the boundary S matrix [32] does depend on 
the boundary parameter. Therefore, the corresponding NLIE may have a Pj,(^) term which 
also depends on the boundary parameters. 

Another interesting problem is to identify the specific feature of the spin-1 models which 
leads, in the continuum limit, to supersymmetry. It may be possible to construct an inte- 
grable open spin-1 chain whose continuum limit is integrable but not supersymmetric, such 
as models considered in [29] . We hope to be able to address some of these questions in the 
future. 
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A Range of boundary parameters 

We derive here the restrictions on the ranges of the boundary parameters given in the main 
text (2.10), (3.9). 

A.l Spin-i 

Our argument depends on two assumptions. First, the imaginary part of Ina(x) is monoton- 
ically increasing. Second, there should be no holes on the real axis (except for a single hole 
at the origin) . The first assumption suffers from exceptions for some values of 7 if the system 
size is small. However, we regard this as a finite size effect and expect our assumption to be 
valid in the thermodynamic limit. 

We consider first the repulsive regime < 7 < |. Let < a < |. We then define two 
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functions, 

n I X 1, / sm\i{x-ia) 1 sinh^x -la) 

^ ' I ^ s\-nh{x + ia)' +v w ^ ^ s\-nh{x + la) ' 

where we adopt the convention for the logarithm such that its imaginary part is restricted 
to [— 7r,7r). These functions are essentially the same but for the choice of phase (Figure 5). 
We choose the branch of (a;, a) such that a branch cut line emerges from ia and goes to 




Figure 5: The functions (left) and (right). 



positive infinity, and another line starts from —ia and goes to negative infinity, as shown in 
Figure 6. 











i-ia 


< 



Figure 6: The branch cut lines for 

In the rest of the range, | < a < tt, we define 

'i?_(a;, a) = — 'j?_(x, TT — a) , 'i?+(a;, a) = — 't?+(x, tt — a) . (A.l) 

Recalling (2.5), we consider 

-Loga(a;) = V d_{x, ^) + d+{2x, 7) + iV V i?_(a; - aA, J) 

M 

-Y,[d^{x-Vk,i)+^-{x + Vk,i)] , (A.2) 
fc=i 
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where M denotes the number of BAE roots. By Log, we mean the logarithm whose imaginary 
part is not restricted to [— 7r,7r). Strictly speaking, a is defined in the upper half-plane. We 
assume that it is also well-defined on the real axis. Note that |Loga(0'*') = '»?+(0"^, 7) = — tt 
for < 7 < 7r/2. 

With the help of the identities 

TT 

'i?_(oo, a) — TT — 2a , '*^+(oo, a) = —2a , < a < — , 
one easily derives the asymptotic value 

— Loga(oo) = 2(A^ - M + 1) - ^{2N -AM + 2 + H+ + H) (A.3) 

in TT 

for H± in the in the interval (2.9). We denote the integer part of the above as n, 

n = 2{N-M+l)-6, 

where 

5 = [-(2A^ -AM + 2 + H+ + H_)\ + 1 , 

TT 

and [x\ specifies the integer part of x. When a; is a root or a hole, 4j:Loga(a;) is an odd 
integer. The range of ^^Loga is set so as to be able to accommodate M roots. As already 
noted, there should be no holes (except for a single hole at the origin). Thanks to the 
assumption of monotonicity, these conditions determine the range of n , 

2M-l<n<2M+l. (A.4) 

This is equivalent to having S = 2,3. The ground state selects M = N/2, and this leads to 
the conclusion (2.10), 

-r-r STT 

- <H+ + H_ + 2<—. 
7 7 

In the attractive regime 7r/2 < 7 < vr, some functions in (A. 2) must be rewritten in terms 
of 7 = TT — 7 which is now in the range (0,7r/2). Thanks to (A.l), this only results in an 
extra 27r on the RHS of (A.3). In this case, however, 4Loga(0+) = tt, which cancels this 
extra 27r in the argument for the possible values of 5. We then reach the same conclusion on 
the range of i?±. 

A.2 Spin-1 

The spin-1 case suffers from an extra complication due to the absence of an auxiliary function 
which precisely encodes the "branch cut integers" such as a and a in the spin-1/2 case. 
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Previous studies [25, 28] nevertheless give this interpretation to lab{x). It roughly encodes 
the information of the two-string center. We assume this here, and utilize b in place of a in 
the above argument. 

We consider first the repulsive regime, < 7 < |. The following representation is 
convenient for our purpose (recall Eqs. (3.1), (3.13)), 

^mLogb{x) ^ Nj2['^-i^ - ?) aA, ^)] 



- J2 [^-(^, Va-l)+Mx,rj^ + ^)]+ ^+(2x, 27) 

M S3 



k=l 

+ 9mln(l + ai(x-^)), (A.5) 

where we define a\{x) = li{x)/l2{x). In the above, i?+(2a;, 27) should be understood as 
— '(9+(2a;, tt — 27) if 7r/4 < 7 < n/3. The above expression is valid for '^mx = +e. 

We assume that the ground state is given by a sea of 2-strings with slight deviations, 
\'^mvk\ =7/2 + 6, and M = N. Then a careful analysis concludes 

Z7 

ln(l + ai{x + ie- —))\^^^ ^ ^ - [r]+ + rj-) + Stt , 

where S is an integer satisfying \ j—{rj++r)-)+STr\ < it. This leads to the following expression 
of the asymptotic value. 

Log 6(00) = 3(77+ + rj- — ^) — An + 2nN + 5n , 

where, for simplicity, we have further restricted the boundary parameters to the interval 
7r/2 < 7y± < TT - 7/2. 

For the spin-1 case, we must locate N/2 2-string centers on the real positive x-axis. Then 
an argument similar to the one for the spin- 1/2 case leads to 

5 Stt 

TT < 77+ + 77_ - 7 + -TT < y. 

One can then easily draw a conclusion that 5 = 1 is the only consistent choice, and that the 
constraint (3.9) must be imposed, 

27r 47r 

y + 7< 77++ 77_<y+7. 

A little modification leads to the same conclusion for the attractive regime 7r/3 < 7 < 7r/2. 
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B Boundary terms 

We provide here details in the computation of various boundary terms appearing in the 
NLIEs. 



B.l Spin-|: Crik) 

In the spin- 1 calculation, the quantity Crik) is defined by (2.25) 

CT{k) = G{k) C{k) + D{k) . 



(B.l) 



Since B'^^^x) and 0(a; ± ^) are ANZ near the real axis, the contour integral on Ci in (2.18) 
can be changed into — C2, and C{k) becomes 



C{k) 
5{k) 



-D{k) - 5{k) , 



I [ln/x(x)]"e^*^^ = -ik I [ln/x(x)]'e''=^ = 27rA; . 
Jc Jc 



It follows that 



CT{k)^D{k) 1-G{k) -27rkG{k) 
We now evaluate D{k), which is defined by (2.21), i.e.. 



(B.2) 
(B.3) 

(B.4) 



D{k) 



sinh(a; + ^) sinh(a; + ^) sinh^(x - A + f ) sinh^(x + A + |) 



sinh(a; - ^) sinh(a; - ^) sinh^' (^ - A - ^) sinh^' (x + A - ^) 



X 



sinh(2a; + 27) 



sinh(2a; — ^7) 

Making use of the identities (2.26) and (2.27), we obtain 



2w 



iAk\ 



g(2 ^)'= _ e"^*^ 



( X — 2L]U —Ik- 



(B.5) 



We have assumed here that the boundary parameters Hj- are in the domain (2.10). 
Using these results we can simplify the first term in Eq. (B.4). From (2.24), 



7 

ke 2 



4cosh(f )sinh((7-7r)|) 



ke^ cosh ^ 



2cosh(f )sinh((7-7r)|) ' 
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Combining these, 
1 



2'Kk 



D{k)[l-G{k)\ = 



iVcos(AA;) cosh(^)sinh((7- f)|) 



cosh(f) cosh(f )sinh((7-7r)|) 

sinh {{^H+ - 7^)1) + sinh ((7-^- - 7r)|) 
2cosh(f )sinh((7-7r)|) 



Finally, taking into account also the second term in (B.4), we conclude that Cr(A;) is 
given by 



CT{k) 



^^^J n cosjAk) sinh {(-fH+ - 7r)|) + sinh ((7//, - 7r)|) 
I cosh(f ) ^ 2cosh(f)sinh((7-7r)|) 



cosh(f)sinh((27-7r)|) 
cosh(^)sinh((7-7r)|) 

B.2 Spin-1: C(A;) 

In the spin-1 calculation, the quantity C{k) is defined by (3.36) 

C{k) = -G{k) Dg{k) + G2{k) LT'ik) + LC^{k) - e-'rl^'ik) + -J^ 



(B.6) 



'yk 

e '2 + e 2 



(B.7) 



e"^ + e"^^) L(P"{k) - Ln"{k) , (B.8) 



The quantity Dq{k) (3.30) can be evaluated as follows: 

Lifl{k) = 2T:^2{k)e-^^ + (y + e"^) LBM" {k) + 
-£i^(A;) = 27rV'2(A;)e('^-t)'=+ (^e^+e-^)LSH"(A;)+ [e(^-'^)'= + e(2-'^)*^] L0^'(A;) + £i?'(A;) , 
which gives 



Prom the definition of f{x) (3.5), we can derive 

L?(fc) = 



2 -)- e 2—6^2 'rj«_gV2 '^J 



[ln/i(a;)] e 



C 



/ dx< 




IC2 





sinh(2a; - 2^7) sinh(2a; + 2i7)5(+)(a; - —)B^'\x + — ) 



Akx 
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It follows from the identities (2.26) and (2.27) that Lf"{k) is given by 

L^'ik) = 27rV^2(A;)[e-^'= + e(^-5)'=]+ [e-^L^W''(A;) + e^LB(-)"(A;) 

Similarly, from the definition of Cb{x) (3.15), 

Lq;{k) = -2'KiP2{k)e^^-^^' + e'^LB(+y'{k)-[e^ + e-^jLB(-y'{k) 

+ 27r^(A;)iV(e^^'= + e"'^") [e"^ - e(i-")*^ - e(^-")] + e-^5?(/t) . 
Substituting these results into (B.7), we obtain 



(B.9) 



C{k) 



2 cosh f 



LBM" (k) - LB(^" (k) 


/ irk irk \ 

( e 2 — e 2 1 




e(f-^)^ - e(^-5)'" 


1 e 2 + e 2 1 



2Tiil)2{k) ( 1 — e '2'° I I — e '2 



( e 2 + e 2 1 



;f-7)fe_p(7-f)fe 



+ 27rfc(G(fc)+ ^ ^ ^ ^ , 
e 2 + e 2 



where we have also used the result (B.3). From the definition of B'^^\x) (3.3), 

LBM"{k) = 27rV^(fc) [e(''+-^)'= + e^"— "^'^j , 
LB^"{k) = 27rV'(A;) (e-^+^^ + e-^-*^) . 

We conclude that C{k) is given by 

[smli((//+-f)/,-)+smli((//_-f)fc)] 



(B.IO) 



C{k) 



2Tik< N 



2 cosh ^ 



+ 



2coshf sinh((f -7)A;) 



+ 



cosh^sinh((37-7r)|) 
coshf sinh ((27-7r)|) 

B.3 Spin-1: Cyik) 

The quantity Cy{k) is defined by (3.42) 

Drik) 



(B.ll) 



e 2 + e 2 



:b.i2) 
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Prom the definition of DT{k) (3.40) and the results (B.8), 

= 27rV'2(A:)[e-^ + e^^-t)'^] + (l + g-^^^) LBW"{k) + (l + e^^^) LB^"{k) 

+ e'^C4Ji!'{k)-e-'^Liji/'{k). 
It follows from (B.12) that Cy{k) is given by 



Cy{k) = 277^^2 (A;) 



e 2 + e*- 2 2 
+ e 2" 



e-'^LBi+)"(k) + e'^LB(-)"(k) 



iAk I iAfc\ 



+ 2TTip{k)N{e''''' + e 



2 ^ + e 2 



e 2 + e 2 



+ 27rij2{k)-Lf"{k) 

= 47rA;G'2(-A;) , (B.13) 
where, in passing to the last line, we have used the results (3.35), (B.3) and (B.9). 

C Integration constants 

We explain here how to determine the integration constants in the NLIEs. 



C.l Spin-i 



The main idea is to carefully consider the limit x — > oo. In this limit, the NLIE (2.29) 
becomes 



In a{oo) = G{oo) [ln74(oo) — ln74(oo)] + i Phdryi^o) + iirC , 



(C.l) 



where C is the constant which is to be determined. The contribution from the driving term 
is a multiple of 27ri (since N is even), and can therefore be dropped. Prom the definition of 
a{x) (2.5) and the fact that M — N/2 for the ground state, we readily obtain 

a(oo) = e''^ , u = -f{H+ + H_ + 2). (C.2) 

We define the branch of the logarithm such that \Qm lnx\ < tt. It follows that 

Ina(oo) = i {u — 27rm) , (C.3) 
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where m is an integer such that 

-n < (J - 27rm< n . (C.4) 

Moreover, 

A(oo) = 1 + a(oo) = 1 + e''^ = 2e^ cos ^ = 2e'(t+^'^) I cos ^ | , (C.5) 



where S is given by 



S={ ° ■ . (C.6) 

' ±1 : COS I < ^ ^ 



Hence, 

In^(oo) = In 2| cos^l + i + 7r5 - 27™) , (C.7) 

and ln74(oo) is obtained by complex conjugation. Here n is another integer such that 

-TT < ^ + 7r5 - 27m < TT . (C.8) 

We now substitute into (C.l) the expressions for Ina(oo) (C.3) and liaA{oo), ln/l((X)) (C.7), 
as well as the results 

= 2V-7) ' P>"^ry{oo) = i^H^ + + 47 - 47r) , (C.9) 

which follow from (2.30) and (2.31), respectively. Solving for a;, we obtain 

uj = + 7if_ + 27r (C + 5 + 2m - 2n - 2) + 27 (-C - 25 - 2m + 4n + 2) . (C.IO) 

Comparing this result with the definition of a; in (C.2), and assuming that C is independent 
of 7, we obtain a pair of equations, 

C + (5 + 2m-2n-2 = 0, -C-25-2m + 4n + 2 = l, (C.ll) 

which imply 

5 = 2n-l (C.12) 

and 

C = 3-2m. (C.13) 

The relations (C.12) and (C.8) imply that < u < An. In fact, since 6 can be only or 
±1, (C.12) implies 6 = ±1. It follows from (C.6) that u is further restricted to the interval 

7r<a;<37r. (C.14) 

Finally, (C.4) then imphes m — 1, which determines C through (C.13), 

C=l. (C.15) 

Note that the definition of uo (C.2) together with (C.14) imply the domain of boundary 
parameters quoted in the text (2.10). 
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C.2 Spin-1 



In the limit x — > oo, the spin-1 NLIE becomes 

ln6(oo) = G{oo)[lnB{oo)-lnB{oo)]+^lnY{oo)+iPbdry{oo)+i7iC, (C.16) 
Iny(oo) = ^ [ln5(oo) +ln5(oo)] +i7rCj/, (C.17) 

where C and Cy are the constants which are to be determined. As in the spin-| case, we 
assume that N is even, and therefore drop the contribution from the driving term. Prom the 
definition of h{x) (3.12) and y{x) (3.21) and the fact that M — N ior the ground state, we 
obtain 

6(00) = e'-'" + e'*" , y (00) = 1 + e^"' + e"'^'" , (C.18) 
where now cu is defined as 

a; = 77+ + 77_ -7. (C.19) 
We rewrite the expression for 6(00) as 

6(00) = 263^"^ cos u; = 2e^'^^+''^'^\cosu}\ , (C.20) 

where Si is given by 

5, = | ° ■ . (C.21) 
I ±1 : cos a; < 

It follows that 

In 6(00) = In [2 cos a; |] +i{3uj + ttSi - 27rm) , (C.22) 
where m is an integer such that 

-n <3u; + nSi - 27rm < tt . (C.23) 

Moreover, 

B{oo) = 1 + 6(00) = 1 + e'^'" + e^^ = e'^il + 2cos2a;) = 2e*(''^+'^^2)|l + 2cos2u;| ,(C.24) 
where 82 is given by 

^ ±1 : l + 2cos2a;<0 ^ ^ 
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Hence, 

InS(oo) =ln|l + 2cos2tu| + i(2a; + 7r52 -27rn) , (C.26) 
and ln5(cxD) is obtained by complex conjugation, where n is an integer such that 

-n <2u; + 77(^2 - 2?™ < tt . (C.27) 

From (C.18) we also see that 

y(oo) = 1 + 2cos2u;, (C.28) 

and therefore 

F(oo) = 1 + 2/(oo) = 2 + e'*" + e-'*" = (e^" + e-^' f = |2 cosa;|' . (C.29) 

We now substitute into the first NLIE (C.16) the above expressions for ln6(oo) (C.22) 
and InS(oo), lni?(oo) (C.26), as well as the result 

G(oo) ^ ^ (C.30) 

which follows from (3.34), and the expression for Pb^ryioo) (5.5). Solving for cu, we obtain 

cu = r/+ + r/_ + TT (C - 5i + ^2 + 2m - 2n) + 7 (-2C + 26^ - 3^2 - 4m + 6n - 3) . (C.31) 

Comparing this result with the definition of cu in (C.19), and assuming that C is independent 
of 7, we obtain a pair of equations, 

C -Si + S2 + 2m-2n^0, -2(:7 + 25i - 3(^2 - 4m + 6n - 3 = -1 . (C.32) 

These imply that S2 = 2n — 2, which is an even number. But since S2 can be only or ±1, 
this imphes 

62^0, n = 1 . (C.33) 

It follows from (C.32) that 

C = 5i - 2m + 2 . (C.34) 

The relations (C.33) and (C.27) imply that f < a; < ^. Hence, 

5i = ±1 . (C.35) 
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It follows from (C.25) and ^2 = that uj is further restricted to the interval 



27r 47r 



(C.36) 




Finally, (C.23) then implies 6i — 2m 



3, which determines C through (C.34), 



C= -1. 



(C.37) 



Similarly, substituting into the second NLIE (C.17) the results (C.26), (C.28), and remem- 
bering that 62 = 0, we immediately see that 



Note that the definition of uj (C.19) together with (C.36) imply the domain of boundary 
parameters quoted in the text (3.9). 
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